Major structures of the inner Galaxy delineated by 6.7-GHz methanol masers 
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ABSTRACT 

We explore the longitude- velocity distribution of 6.7-GHz methanol masers in the 
context of the inner structure of our Galaxy. We analyse the correlation in velocities 
within this distribution and identify density enhancements indicating large-scale regions 
of enhanced star formation. These are interpreted as the starting points of the spiral 
arms and the interaction of the Galactic bar with the 3-kpc arms. The methanol 
masers support the presence of a long thin bar with a 45° orientation. Signatures of 
the full 3-kpc arm structure are seen, including a prominent tangent at approximately 
—22° Galactic longitude. We compare this distribution with existing models of the gas 
dynamics of our Galaxy. The 3-kpc arm structure appears likely to correspond to the 
radius of corotation resonance of the bar, with the bar on its inner surface and the 
starting points of the spiral arms on its outer surface. 

Subject headings: masers — physical data and processes, stars: formation — stars, 
Galaxy: structure — The Galaxy 



1. Introduction 



Studying the structure and dynamics of the Milky Way is complicated by our location within 
it, and the resultant obscuration along the line-of-sight, but nevertheless it has established that our 
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Galaxy has a grand spira l structure, with spiral ar ms, interior structures such as the 3-kpc arms 
( van Woerden et alj 1957 ; Dame h Thaddeus 2008 ) and one o r two Galactic bars ( Blitz k, Spergel 



1991 



Gerhard! 120021 . references therein; iBenjamin et al.ll2005l ). Early descriptions of the structure 



of our Galaxy were inferred from star counts and positioning of objects with photometric distances, 
but current revisions have been largely based upon the kinematics of atomic and molecular tracers, 
such as Hi and CO. The sinusoidal distribution of velocities with Galactic longitudes has established 
the approximately circular and symmet ric rotation of the Galaxy and both tracers have been 
used to constrain dynamical models (e.g. 



Binnev et all Il99ll : IWeiner fc Sellwoodl 1.1999 



Peters 



1975 



Fux 



Burton fe Liszt 1978; Liszt fc Burton 1978 



1999). However, there are limitations to this 



approach. Hi is detected throughout the Galaxy in such great abundance that only the highest or 
lowest velocities along a given line of sight (the terminal velocities and spiral arm tangent deviations) 
can be used to constrain dynamical modelling. CO is less diffuse and more concentrated in the 
major structural features of our Galaxy, but surveys of CO are limited by a lack of sensitivity to 
material on the far side of the Galaxy. An alternative tracer, the population of 6.7-GHz methanol 
masers, has the advantage of only being observed towards regions of high-mass star formation 
(|Pestalozzi et alJbood : llVlinier et all 1200.4 ISridharan et alJbood : [Walsh et al.ll20o4 IXu et al.ll2008h 
and the strong and narrow emission of these masers is well correlated with their systemic velocity , 



the central maser velocity typically within 3-5 km s 1 ( Szvmczak et al 



They are also found in abundance throughout the Galaxy (IPe stalozz i et al 



200 



7 



Pandian et al.1 120091 ) 



20051 ) and have been 



detected to Galactocentric distances of 13.5 kpc (jHonma et al.1 120071 ) . Inner Galaxy structures 
have been sh own to be traced by 6.7-GHz methanol masers: 45 sources are associated with the 
3-kpc arms (IGreen et al.l 120101); 7 sources with the h igh velocities associated with Galactic bars 



(lCaswelllll997l ; ICaswell et al.ll2010l : iGreen et al.ll2010l ); and 11 with Sagittarius B2 ([Caswell et al 



2010 . references therein). The 6.7-GHz methanol masers are excellent measures of the kinematic 
behaviour of the structure of our Galaxy within the terminal velocities outlined in Hi and CO. 

In this paper we explore how 6.7-GHz methanol masers trace the major star forming structures 
of the inner Galaxy, specifically: the Galactic bar, the 3 kpc arms, and the origins of the spiral arms. 
All of these structures are believed to be confined within ~4 kpc of the Galactic centre. Assuming 



a solar distance of 8.4 kpc (jGhez et al.l 120081 ). the inner 4 kpc of the Galaxy is geometrically 
contain ed within the Ga lactic longitude range —28° to +28°. The Methanol Multibeam (MMB) 
survey (IGreen et al.ll2009l ). which observed from —174° to +60° with the Parkes Radio Telescope, 
detected ~550 masers within longitudes ±28°. The MMB had a 3 a sensitivity of 0.7 Jy, covered the 
latitude ran g e of ± 2° and had a velocity coverage encompassing the extent of CO emission seen in 
Dame et al.1 (120011 ). The MMB catalogue therefore represents the most complete survey of current 
high-mass star formation traced by 6.7-GHz methanol masers in the region, providing a census 
of the regions which includes their kinematic behaviour. Through examining the distribution and 
densities of these masers in longitude and velocity, in comparison with the spiral arms, we attempt 
to delineate the inner Galactic structure. 



-3 - 



2. 6.7 GHz methanol maser distribution and density enhancements 

The longitude- velocity (l,v) distribution of 6.7-GHz methanol masers between longitudes ±28° 
is shown in Figure[TJ It is immediately apparent by eye that the maser distribution is not entirely 
random, but contains clumps of sources on the sub-degree scale over a few kilometres per second, 
together with structures on the scale of several degrees over tens of kilometres per second. The 
former is due to multiple sources within molecular cloud complexes, whilst the latter is a signature 
of Galactic scale structures. A structure function of the l,v distribution formally demonstrates the 
presence of this structure, quantifying the velocity correlation for given angular separations. We 
define the second order structure function of velocity as a function of angular separation: 



SF vel (r) = ([vel(x) - vel(x + r)] 2 ) 



(1) 



where x is the angular position of a maser site and r is the angular separation between a pair 
of maser sites. We take the ensemble average of measurements with the same range of angular 
separation (denoted by the angular brackets) . The resulting structure function is shown in Figure [2] 
for both the masers and a Monte Carlo simulation together with statistical errors. The Monte Carlo 
simulation is a random distribution of masers in the plane of the Galaxy with longitudes between 
±28° and Galactocentric distances between 3 and 13.5 kpc. The masers were assigned velocities 



based on a Galactocentric solar distance of 8.4 kpc (jGhez et al 



2008; 



Gillessen et al 



2009) and flat 



rotation curve with a circular rotation of the Sun of 246 kms -1 ( Reid et al. 20091 : Bovy et al. 2009 ). 
We generated 482 sources by this method, to match the observed quantity for the l,v domain, i.e. 
not including the 63 masers mentioned in the introduction associated with the 3 kpc arms (45), 
Galactic bar (7), and Sagittarius B2 (11). A power law fit to the structure function gives a gradient 
in the log-log plane of +0.09±0.01 for the masers and — 0.02±0.01 for the Monte Carlo simulation 
(with the errors the formal errors in the fit). The consistent positive gradient of the structure 
function of the data is indicative of a dumpiness of structure on large scales (both angular and in 
velocity) and a correlation on small scales. In comparison, the approximately flat structure function 
of the Monte Carlo simulation demonstrates the same level of velocity correlation at all scales, i.e. 
no scale-specific structure. The measured structure function departs from the simulation, within 
the errors (both formal and statistical), at scales less than 0.03° (with correlated velocities up to 
30 kms -1 ) and at scales greater than 3° (with correlated velocities up to 50 kms -1 ). Smoothing 
the distribution of sources on the small scales (by counting all sources within 0.03° and 30 kms -1 
singularly) and binning the data on the large scales we produce the l,v density distribution seen in 
Figure [3l The non-zero mean of the bins is 3.7 sources, giving a Poisson noise of 1.9 sources. At 
a 5cr confidence level, statistically significant density fluctuations are therefore those with 9.7 or 
more sources, and these are highlighted in Figure The density enhancements are also apparent 
in the distribution of number counts with longitude, shown in Figure [H 



The majority of the Galactic scale structure seen in the 6.7-GHz methanol maser l,v domain 
is interpretable as the spiral arms. Figure!]] overlays the masers on the loci of the spiral arms 
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comm only adopted in analyses of struc ture in the l,v domain: the fou r arm shape of lTavlor &c Cordes 
(|1993l ). effectively the same as that of iGeorgelin Georgelinl ()1976l ). transferred to the l,v domain 
using the rotation curve of brand k Blitzj (|l993l ~ We adopt the solar rotation and radius as 
mentioned in the previous paragraph and use these throughout the discussion. The choice of 
rotation curve, solar parameters and spiral arm shape have minimal effect, but these are addressed 
in detail in the appendix. Sources already associated with star forming regions towards the Galactic 
centre (Sgr B2, the Galactic Centre Zone and sites within 3.5-kpc of the Galactic centre, see 



Caswell et al.l (|2010l ) for details) and the parallel sections of the 3-kpc arms between longitudes ±15° 



are distinguished by different symbols. With these excluded from consideration, then incorporating 
an arm thickness of 1 kpc with a velocity tolerance of ±7kms _1 , we find the spiral arms account 
for 79% of the masers. The Monte Carlo simulation (from 100 realisations) gives an average (both 
mean and median) number of associated sources of 69% with a standard deviation of 2%, implying 
the association of real masers is statistically significant (none of the Monte Carlo simulations exceed 
72% association). We chose an arm thickness of 1 kpc based on estimates of the inter-arm separation 



and typical widths of the arms in models (see fo r example I Gomez &i Cox 



2004; 



Sewilo et al 



2004; 



McClure-Griffiths et al.ll2004l : iLevine et al.1 120061 ) . Varying the arm thickness between 0.5 and 1.5 
kpc results in the level of association varying between 49 and 78%. Our velocity tolerance of 
±7kms~ 1 is chosen based on the kin ematic uncertainty of the masers in re lation to the high-mass 
star forming region they are tracing (jReid et al.ll2009l ; iPandian et al.ll2009i ). Varying the velocity 
tolerance between 5kms _1 and lOkrns -1 results in an association of between 58 and 80%. The 
unassociated sources (approximately 20%) lie almost exclusively in the longitude regions 18° to 28° 
with velocities >100kms _1 and between —19° and —24° with velocities between and — 130 kms -1 . 
These sources will be discussed in the following sections. 



3. Structural interpretation of l,v distribution 

The /^distribution raises the question: what are the causes of the statistically significant over- 
densities — are they related to Galactic structure or random fluctuations in star formation? Regions 
of star formation are contained within Giant Molecular Clouds (GMCs) which have physical scales 
of 10s of parsecs. With the high-mass star formation traced by the 6.7-GHz masers, we are looking 
at regions at typical distances of >2 kpc (the l,v distribution within ±28° suggesting mostly >4 
kpc). At these distances most GMCs subtend angular scales smaller than 1° . The coherent dense 
regions of sources seen in Figure [H and listed in Table [TJ are on scales of >1° (highlighted by their 
presence post-smoothing in Figure [3]). This makes random fluctuations, unrelated to Galactic scale 
processes, an unlikely cause of the increased star formation. This section analyses the 6.7-GHz 
methanol maser l,v distribution in relation to the major structural features of the inner Galaxy. 
We begin with the Galactic bar, then the 3-kpc arms, and follow with the starting points of the 
spiral arms, or the spiral arm 'origins'. Excluded from the discussion are masers at high latitudes, 
as these are unlikely to lie w ithin 4 kpc of the Galactic centre. With a scale height of star formation 



estimated to be ~30— 40 pc (jFish et al.ll2003l ). any masers whose latitudes imply a Galactic height, 
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Fig. 1. — Longitude- velocity 'crayon' plot showing the distributi on of 6.7-GHz m e thano l masers within 
±28° overlaid with e xample spiral arm loc i (the spatial pattern of iTavlor fc Cordesl (|l993l) . transferred to 
the l,v domain by the Brand fc Blitz! (1993) rotation curve). Coloured loci are the spiral arms defined by the 
model, grey shading incorporates an arm thickness of 1 kpc and a velocity tolerance of 7kms _1 . Yellow loci 
represent the Perseus spiral arm; Purple - Carina-Sagittarius; Orange - Crux-Scutum; Green - Norma. The 



blue lines delineate the region identified in CO emission as the 3-kpc arms bv lDame fc Thaddeud (|2008l ). 
Crosses s how 6.7-GHz methanol masers of the MMB survey. Circ les are masers associated with the 3-kpc 
arms (see Green et al. 20091 : Caswell et al. 2010 ; Green et al. 2O10l for details). Diamonds are masers which 



are in terior to the 3-kpc arms, primarily candidates for belonging to the Galactic Centre Zone (jCaswell et al 
20101 ). Triangles are masers associated with the Sagittarius B2 complex. Stars are masers associated with 



the Galactic bar. Crosses enclosed in squares are masers with high latitudes (and therefore likely to be 
closer to us than 4.5 kp c, see Section [2]). The starting points of the spiral arms have been adjusted from 
Tavlor fc Cordesl (|1993l ) to match the discussion of Section 13.31 The red brackets highlight the longitude 
range of the dense ridge of sources discussed in the text. 
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Fig. 2. — Velocity structure function showing the correlation in velocity for given angular separations of 
maser sites. The data is binned in equal logarithmic bins (0.3) and compares all the MMB data between ±28° 
longitude (red) with a Monte Carlo simulation (black) . The Monte Carlo simulation has been systematically 
offset in angular separation by +0.05 for clarity of comparison. The error bars represent the statistical 
(Poisson) error in the mean of each bin. 
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Fig. 3. — Longitude-velocity 'blockbusters' plot showing the density of masers within the l,v domain. As 
described in the text the distribution of sources seen in Figure [1] has been smoothed on scales <0.03° and 
then binned on scales >3°. Hexagonal bins were used as they reduce the visual distortion in the distribution 
caused by binning. Fhe colour-scale represents the number of sources. Black lines delineate statistically 
significant dense regions, at the 5 (J po isson level. Labels correspond to enhanced star formation at the ends 
of the Galactic bar and the starts of the spiral arms. 
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Fig. 4. — Distribution of sources with longitude in 2° bins. Median number count per bin is 18 (red line). 
Errors represent Poisson statistical error margins. Beyond the limits of the figure, existing maser data suggest 
the number counts decline to (or below) the average beyond —28°, but show a peak at positive longitudes, 
beyond +28°. 
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Table 1: Parameters of density enhancements of 6.7-GHz methanol masers in the longitude- velocity 
domain. For explanation of source density see main text. Densities with + are the sum of more 
than one enhancement, confused by multiple passes of spiral arms along the line of sight. The last 
column lists Galactic structures which are believed to be responsible for the enhancement. 
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z, greater than ±60— 80 pc at a heliocentric distance of 4.4 kpc can be excluded from inner Galaxy 
considerations. This corresponds to latitudes 0.8° to 1.0° and we therefore adopt 0.9° as a cut-off. 
Within the longitude region ±28° there are 29 masers with latitudes in excess of 0.9° from the 
plane. 



3.1. The Galactic bar 



The Galactic bar has usually been regarded as either a triaxial boxy bulge (or 'peanut'), princi- 
pall y on the basis of Cosmi c Background Explorer (COBE) infrared observations, or as a long thin 
bar. iBlitz Sz Spergell (|199ll ) speculated that the Gal axy might contain two bar components and re- 



cent consensus is that both these types of bars exist (jHammerslev et al. 



2000 



Benjamin et al. 



2005 



Cabrera-Lavers et al.ll2008l ). The orientation of the bar with respect to the Sun-Galactic centre line- 
of-sight is thought to lie between 14° and 45°, but many of the lower inclination angles are for boxy 
bulge bar models, whilst the larger angles are for long thin bars (and confusion has been added with 
the in terch angeable and sometimes ambig uous use of the bulge and bar terms, see discussions oflNg 
(|1998l ) and lLopez-Corredoira et al.1 ([1999)). This has led to a likely picture of a triaxial bulge bar 
inclined at 20—30° and a long thin bar at ~45°. This combination is shown relative to the 3-kpc 
arms in Figure [5j Primarily on the basis of models based on infrared (IR) data, the radius of corota 



tion r esonance of the bar (of the order of 1.2 times the se mi-major axis of the bar (jElmegreen et al 



1996; lEndmaier & Gerhard 12009 : iButa k Zhand 120091 )) is estimated to be 4±0.5kpc (iGerhard 



20021 . references therein). In galaxies t his radius marks the transition from ellip tical orbits follow- 



ing the bar to circular orbits beyond (jContopoulos Papavannopoulod Il980l ). As noted in the 



introduction, several 6.7-GHz methanol masers have been associated with the Galactic bar. Where 



formation is expected (e .g. iFux 



1999 



Englmaier & Gerhard 


1999: 


Lopez-Corredoira et al. 


2001 



Binney &: Tremaindl2008l ). This means that a density enhancement in either the near or far 3-kpc 
arm 6.7-GHz methanol maser population could indicate the interaction w ith the end of the bar. 
Such an enhancement has been suggested in the far 3-kpc arm population (IGreen et al.ll2010l ) and 
we now address both this and possible near-side interactions. 



3.1.1. Density enhancement at the far end of the bar 



The 6.7-GHz methanol maser population tracing the 3-kpc arms exhibits a high density of 
sources within the far 3-kpc arm between longitudes —9° to —15° (see Figure[3]and Table[T]). Within 
this dense region, 8 sources with velocities >0kms _1 partially overlap with the extrapolated outer 
Galaxy components of the Carina-Sagittarius and Crux-Scutum arms (at distances > 15 kpc), and 
6 sources with velocities between and — 30kms _1 overlap with the Carin a-Sagittariu s arm (which 
passes between the Galactic centre and the Sun at a distance of <1.1 kpc lSato et al.1 (120101 )) . 
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The distances of these arms will have an impact on the flux density distribution of the sources, 
an effect that can b e estimated if we assume a simple power law luminosity distrib ution (e.g. 



Pestalozzi et alj|2007l ) above the sensitivity limit of the MMB (0.7 Jv. iGreen et al.ll2009l ): at 15kpc 



heliocentric distance only the brightest sources would be detected, whilst at 1 kpc heliocentric 
distance very faint sources would be. If the far 3-kpc arm sample of sources was significantly 
biased by distant sources, we would expect a shallow distribution of peak flux densities, whilst if it 
was biased by nearby sources, we would expect an asymmetric distribution with a peak at very low 
peak flux densities. We do not see either, instead finding an approximately Gaussian distribution 
with a median of 2.4 Jy, implying a population dominated by intermediate distances. Furthermore, 
the far spiral arms are not well traced by Hi or CO at these longitudes, so their location is less 
precise. Regions of high-mass star formation are comparatively rare, so without either of the spiral 
arms being tangential, we would not expect there to be a significant populat ion present in either . 
The se factors combine to give us confidence that the dense region identified by ICaswell et al.l (|2010i ) 
and lGreen et al.l (|2010l ) is a product of far 3-kpc arm sources interacting with the bar. 



3.1.2. Density enhancement at the near end of the bar: implications for length and orientation 



Expanding on t he analysis of IGreen et al.l (|2010l ). we assume a bar orientation of 45° (e.g. 
Benjamin et al.ll2005l ). with the semi-major axis left as a free parameter. The density enhancement 
at the far end of the bar (between longitudes —9° and —15°) thus implies a semi-major axis of the 
bar between 2.2 and 4.3 kpc. This indicates that the near end of the bar lies between longitudes 
13° and 30°. We see from Tableland Figure [3] that there are three higher density regions of masers 
within this range: between 10° cind. 14° (1 bin above 5 o~p isson 

); between 18° and 22° longitude (1 
bin above 5 a po i sson ); and between 24° and 28° longitude (1 bin above 5a po i sson ); these correspond 
to semi-major axes of 1.8-2.4 kpc, 2.9-3.4 kpc and 3.7-4.1 kpc respectively. In dynamical simula- 
tions, bars with short semi-major axes (or the assumption of such a mass distribution) are almost 



exclusively associated with small bar orientation angles (e.g. iBinney et al. 



1991 



Freudenreich 



1998 



Babusiaux Gilmorell2005l ) and a short bar with an acute orient ation angle would produce densi- 
ties in the far 3-kpc arm at longitudes smaller than is observed (IGreen et al.1 12010 ). Additionally, 
a shorter bar of this type is associated with the Galactic bulge, but the 6.7- GHz methanol mase r 
population has a narrow latitude distribution unassociated with the bulge (ICaswell et al.ll2010l ). 
Hence the semi-major axis is unlikely to be in the range 1.8-2.4 kpc. The estimate of a radius of 
corotation resonance at ~4kpc implies that a semi-major axis of 3.7-4.1 kpc is also unlikely. 

If, on the oth er hand, we assume the semi-major axis of the bar is fixed at 3.5 kpc (the estimate 
of iGerhardl 120021 and references therein, scaled to 8.4 kpc) and the bar orientation is instead left 
as a free parameter, the density enhancement at the far end of the bar implies a bar orientation 
between 35° and 53°. This range of orientation angles would locate the near end of the bar between 
longitudes 20° and 24°, which overlaps with the high density of masers seen at velocities close to 
+50kms _1 (1 bin above 5a po i sson in Figure[3]). 
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In summary we believe the maser population traces the influence of a long thin bar with a 
semi-major axis of ~3.4 kpc and an orientation of ~45°. This implies that, if both a long and 
short bar exist within our Galaxy, it is the long component which is primarily traceable by (high- 
mass) star formation and is younger. The short boxy/bulge bar was primarily identified in the 
IR observations of evolved stars and lacks 6.7-GHz methanol maser emission, thereby indicating 
it traces an older, more ev olved population. Hence, our observations agree with the implied age 



dichotomy in the bars (e.g. lLopez-Corredoira et al.ll200ll ). 



3.2. The 3 kpc arms 



The (near) 3-kpc arm was originally discovered by Ivan Woerden et al.l (|1957l ) as an absorption 
feature with a ra dial velocity of approximately — 50kms _1 at 0° longitude. It was named by 
Port et al.l (|1958l ) based on the perceived Hi tangent point at longitude —22° corresponding to a 
Galactocentric radius of 3 kpc (when the solar distance is 8.2 kpc). At this longitude the distinct 
feature blends with spiral arm emission. Subsequently the negative longitude tangent has also been 
inferred as the c ause for peaks in the lo ngitudinal b rightness profi l e of radio c ontinuum at low 
frequencies (e.g. iBeuermann et al.1 (119851) a nalysis of lHaslam et al.1 (I198ll. I1982T) 408 MH z data), 



2.4 iim. emission (e.g . Havak awa et al 



1981 



OH/IR star kin ematics (e.e. ISevenst er 1999) and IR 
star counts (e.g. IChurchwell et al.ll2009l ). However the work of lHayakawa et al.1 ()198ll ) showed that 
the peak was not evide nt in higher frequ ency radio cont inuum emission, such as at 1.4 GHz (e.g. 
Mathewson et al.lll962l) and 5 GHz (e.g. lLockmanlll979l ). A positive longitude tangent has never 
been clearly observed. iBanial ()1980l ) inferred a tangent a t 23.6° (with a velocity of ~110kms~ 1 ) 
by assuming a ring-like structure, based on the ring fit of I Cohen &: Daviesl (|1976l ) for the Hi data 
across the longitude range 355° to 6°. They found CO emission at <110kms -1 , but it was not 
distinct from the rest of the Galactic emission. The question remains as to why we do not see 
a clear positive longitude tangent. One possible explanation is that it could be due to greater 
obscuration along the line of sight, with a larger length of the Carina-Sagittarius and Crux-Scutum 
arms present, compared with the negative tangent (where the spiral arms have an orientation almost 
perpendicular to the line of sight). 



Recently iDame k, Thaddeua (|2008l ) discovered the long speculated far 3-kpc arm counterpart 
in CO emission, tracing both arms between ±15 ° longitude. Forty-five 6.7-GHz meth anol masers 



were associated with these parallel l,v structures (IGreen et al 



2009 



Caswell et al.ll2010l ) and one of 



these, the brightest kno wn methanol maser 9.621+0.196, has a distance of 5.2±0.6kpc determined 
by astrometric parallax (jSanna et al.ll2009l ). This places the maser at a Galactocentric distance of 
3.4 kpc, concurring with our expectations (see previous section) for the Galactocentric radius of the 
3-kpc arms. The astrometric observations also showed the source to have a velocity component in 
a radial direction from the Galactic centre (41 kms -1 ) and a velocity component counter to the 
authors chosen model of Galactic rotation (—60 kms -1 relative to a flat rotation curve with circular 
velocity of 254 kms -1 ). 
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3.2.1. A 3-kpc arm tangent 



The methanol maser distribution demonstrates a large region with a high density of sources 
between I = —18° to —28° (Figure [3|). This region comprises a dense ridge (in velocity) of sources 
near I = —22° and a further clump of sources (with a smaller range of velocities) near I = —27° (see 
Figures [T] and H]) . The clump of sources around I = —27° is dominated by the G333 giant molecular 
cloud and is coincident with both the Norma arm (close to its tangent) and the Crux-Scutum arm 
and thus is excluded from the present analysis. The dense ridge near I = —22° however, extends 
over almost 80kms -1 of velocity, most of which does not align with common logarithmic spiral arm 
model loci. The Crux-Scutum arm lies at velocities less negative than — 50kms _1 whilst the Norma 
arm lies at velocities more negative than — 100 km s -1 and less negative than — 50kms _1 , leaving 
masers in the velocity range of —60 to — 85kms _1 unaccounted for. The only arm in conventional 
spiral arm models which could be extended into this region is the Perseus arm (see next section), but 
by itself this cannot account for the high density and range of velocities observed. The implication 
is hence that the ridge feature is a signature of the 3 kpc arm tangent. The prominent ridge 
in the maser distribution also has (les s obvious) c ounte rparts in the l,v d istributions of CO and 
Hi emission (as originally suggested by 



1976 



Port et al 

Downes fc Genzellll980l:ICasweil fc Haynesl 



(|l958T 0. Hn regions (e.g. iGeorgelin fe Georgelin 
19871 ). pre-ultra c ompact Hll region massive young 



stellar objects (e.g. lUrquhart et alJl2007l ) and ammonia emission ((Walsh et al.ll2008l and A. Walsh 



private communication) . 



3.2.2. Fitting an overall structure 



Combining the density enhancement of masers near longitude —22° with the masers already 
assigned to the 3-kpc arms we have a distribution of two parallel structures between ±15°, which 
join at longitudes between —15° and —25°, over a range of ~100kms _1 . Although the ridge 
is a strong signature for a tangent at negative longitudes, the positive longitude tangent is less 
obvious. There are three regions of enhanced densities at positive longitudes (at approximately 
+12°, +20° and +26°, FigureEJ), but they are all confused by spiral arm loci: the region at +12° is 
coincident with loci of every spiral arm; the region at +20° is overlapped by the loci of the Perseus, 
Carina- Sagittarius and Crux-Scutum arms; and the region at +26° is largely coincident with the 
Crux-Scutum arm (and will be discussed in the next section). The region at +20°, which may be 
influenced by the Galactic bar (Section 13. ip . includes a couple of sources outside the spiral arm 
parameter space and potentially extends to a broad velocity range of 15kms _1 to 80kms _1 . 

The two o r iginal interpretations for the 3-kpc arm str ucture were that of an expanding ring 



(jRougoorl 1 1964 : Ivan der Kruitl Il97ll : ICohen fc Davies 



nance feature or elliptical streamlines (jde Vaucouleurs 



1970 



1976 ) or a non-ex panding oval shaped reso- 



Peterdl 19751 ) . Simple examples of these 



two structures are shown in relation to the Galactic bars in Figure [5] and loci are shown in the I, v 
domain for an example of each from the literature in Figure [6l Development of both these interpre- 
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tations was hindered at the time by the lack of a far-side counterpart and the inability to clearly 
trace structure at velocities less than the terminal. The CO detection of lDame Thaddeusl ((2008) 
and the maser distribution allow us to re-examine these structures. The m aser distribution of the 
MMB again suggests an oval structure in the l,v domain (jGreen et al.ll2010l ). Although readily ac- 
counted for by a ring with an apparent expanding velocity, modelling and theoretical interpretation 
favour elliptical orbits/streamlines. Here we investigate this by firstly fitting the 3-kpc arm maser 
distribution with an elliptical model with constant angular momentum and, secondly, a circular 
ring model incorporating a component of radial velocity outwards from the Galactic centre. We fit 
only to the 3-kpc arm maser population as the individual sources towards the tangent points are 
dependent on the spiral arm model and rotation curve (as discussed in the Appendix). 

We fi rst consider the case of an ellipse , with a ratio of semi-minor to semi-major axes of 
0.50-0.80 (|Peterslll975l : ISevenster et al.lll999l ). an orientation between 25° and 50° and a constant 
angular momentum of between 300 and 700kms _1 kpc -1 . We assume the tangential velocity at 
any given point on the ellipse is given by the angular momentum divided by the Galactocentric 
radius of that point. The best fit is found for a semi-major axis of 4.1 kpc, an axis ratio of 0.54 (a 
semi-minor axis of 2.2 kpc), an orientation of 38° and a rotational velocity of 320 km s -1 (a mean 
circular velocity of 106 kms -1 ). This accounts for 96% of the 3-kpc arm sources and is shown 
(with a radial thickness of 0.5 kpc) overlaid on the maser distribution in Figured The ellipse can 
approximate the parallel sections of the 3-kpc arms well, but does not associate the masers at 
positive longitudes (>15°) with high velocities. 

The effect on the level of association with varying parameters of the elliptical ring model 
is shown in the left panels of Figure In this Figure, each pixel within a panel represents an 
ellipse model for a given set of parameters (Semi-major and -minor axes, orientation and angular 
momentum) and the colour scale represents the fraction of masers associated with the model. 
Extending the allowed range of angular momenta does not improve the fit, with the best fits 
confined within a limited range of ~100kms -1 kpc -1 about 300 to 400 kms" 1 kpc -1 , dependent 
on the orientation of the ellipse. Extending the length range of the axes gives much poorer fits 
below 3.5 kpc, but has comparable association at longer radii (although as noted radii larger than 
4.5 kpc are increasingly unrealistic as they result in tangential velocities exceeding the Hi terminal 
velocities). Varying the orientation of the ellipse below 20° produces poorer fits. The effect of 
changing the orientation and the angular momentum is somewhat linked in the l,v domain, with a 
smaller orientation but larger angular momentum producing similar results to a larger orientation 
but smaller angular momentum. Decreasing the ratio of semi-minor to semi-major axes below 
0.5 (i.e. increasing the level of ellipticity) produces a much poorer fit and increasing the ratio of 
axes much beyond 0.8 (making the ring increasingly circular) introduces the need for the radial 
motion of the circular ring scenario. If the unassociated sources described at the start of Section 
[2] are included in the fitting, then the best fit to the combined source sample is for a semi-major 
axis of 4.2 kpc, an axis ratio of 0.53 (a semi-minor axis of 2.2 kpc), an orientation of 28° and a 
rotational velocity of 370 kms -1 (a mean circular velocity of 121 kms" 1 ). This accounts for 91% of 
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the sources. The variation in the fraction of sources associated with the ring models when including 
the unassociated sources is shown in the right panels of Figure [7J 

For the second case of the circular ring, we take radius, radial velocity from the Galactic 
centre, and rotational velocity as free parameters. We assume the ring will start near the end of 
the bar, so allow the radius to range between 2.5 and 4.5 kpc. We vary the Galactocentric radial 
velocity between 40 and 60kms~ 1 in accor d with the range of velocities seen in the CO emission of 



the 3-kpc arms (IDame fc Thaddeusll2008l ). In accordance with the proper motion of 9.621+0.196 
(ISanna et al.l 120090 and our expectation of a lower rotational velocity towards the inner Galaxy 
(e.g. ICaswell et al.ll2010i ). we adopt a circular velocity between 40 and 60kms _1 lower than the 



Galactic rotation (246 km s -1 ). The highest level of source association is found for a radius of 3.4 
kpc, a Galactocentric radial velocity of 48kms~ 1 and a circular rotational velocity of 201 km s -1 . 
This also accounts for 96% of the 3-kpc arm sources and is shown (with a radial thickness of 0.5 
kpc) overlaid on the maser distribution in Figured 

The effect on the level of association with varying parameters of the circular ring model is 
shown in the left panels of Figure (8) As per Figure each pixel within a panel represents a ring 
model for a given set of parameters (Galactocentric radius of the ring, outward Galactocentric radial 
velocity and circular velocity) and the colour scale represents the fraction of masers associated with 
the model. Extending the allowed range of radial velocities does not alter the best fit. Extending 
the range of radii gives a better fit for a larger radius (5.5 kpc) but this is too long physically (HI 
terminal velocities limit the possible length to ~4.5 kpc). Varying the circular rotational velocity 
by a wider range gives marginally better fits (a few %) for higher circular velocities at a given 
radius. Although not formally part of the fit, the ring model also accounts for the tangent sources 
(Section I3.2.ip . the extreme negative sources between -15 and -24° longitude and the extreme 
positive sources between longitudes +15 and +24°. If the unassociated sources described at the 
start of Section [2] are included in the fitting, then the best fit to the combined source sample is for 
a radius of 3.8 kpc, a radial velocity from the Galactic centre of 47kms _1 and a circular velocity 
of 204 kms -1 . However this is a poorer fit to the 3-kpc arm sources between ±15°. The variation 
in the fraction of sources associated with the ring models when including the unassociated sources 
is shown in the right panels of Figure 



3.2.3. Previous models of the structure 



Numerous alternative structure mod els have been previous ly explored. Through an N-body 
model and comparison with OH/IR stars. ISevenster et al.l (|1999l ) proposed an alternative elliptical 
ring. Although this correlated with the near 3-kpc arm known at the time, it is not a good 
match to the far 3-kpc arm or the —22° tangent. This model, which has a constant circular 
velocity rather than constant angular momentum struggles to produce the parallel nature of the 
near and far 3-kpc arms, requiring much longer radii than are reasonable to fit the distribution. 
If, instead of a continuous structure, the 3-kpc arms are considered to be two separate 'lateral' 
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arms, can the maser population be a c count ed f or? Several authors have explored two la teral arms, 
recent examples being lBissantz et al.l (|2003l ) and lRodriguez-Fernandez &: Combes! (|2008l ) (the first a 
smoothed particle hydrodynamics model based on near infrared luminosities of the Diffuse Infrared 
Background Experiment (DIRBE) in comparison with the CO l,v distribution and the second a 
sticky particles model based on star counts from the Two Micron All Sky Survey (2MASS)), but 
neither of these models account for more than about 30% of the 3-kpc arm maser population 
(a third of the level of association of either the expanding ring or the ellipse). More complex 
asymmetric models w hich incorporate severa l arm-like features emanating from the bar fare no 
better: the models of iMulder &: Lieml ((1986) (a two dimensional quasi-st eady state model with 



near a nd far inner Lindblad radius arms and inner ultra-harmonic arms), lEnglmaier fe Gerhard 
(1999) _ a quasi-equilibrium flow model based on COBE/DIRBE near infrared luminosities) and 



Fuxl (|1999l ) (a three dimensional N-body model based on COBE/DIRBE near infrared luminosities 
incorporating the 135 km s^ 1 arm and connecting arms), each associate no more than 30% of the 
3-kpc arm maser s. 

The details of all the models discussed above are summarised in Table [2] and their loci are 
shown in Figures [10] and QTJ The models presented are as published, incorporating the authors' 
chosen best parameters (e.g. pattern speeds, bar orientations, solar distances) and also including 
a lOkms -1 error margin in the loci of the structures. A direct comparison between the types of 
models is complicated by the simple fact that the four arm models cover more of the l,v domain, 
increasing their likelihood of associating sources . How ever, within the elliptical/lateral two arm 
models, that of iRodriguez-Fernandez &: Combed (120081 ) provides the best fit to the sources. This 
model has a large corotation radius of 5—7 kpc, as a result of a slow bar (where that bar is made up 
of both a bulge and longer component). The inclusion of a long bar component, clearly influential 
on high- mass star formation (Section 13.1. 2p . may be responsibl e for this higher associa tion. Within 
the more complex multi-arm models, the resonance features of lMulder &; Lieml ()198 6) provide the 
highest level of association. The lack of significant 6.7-GHz methanol maser emission towards the 
135 kms -1 and 'connecting' arms in these models suggests these may be primarily gas features, not 
undergoing high-mass star formation. This may be due to the broad range of velocities, which both 
features exhibit, indicating highly turbulent material. We similarly do not see 6.7-GHz methanol 
maser emission towards the broad range of velocities of material towards the Galactic Centre (the 
methanol velocity range is ~200kms _1 comp ared with th e CO and Hi spanning ~500kms _1 ) and 
no emission is seen towards Bania's 'Clump 2' (|Banialll977l ). a feature wi th a high range of velocities, 
but it is seen towards the smaller velocity range of Bania's 'Clump 1' (|Green et alJl2010i ). 



3.2.4- Summary of the 3-kpc arm structure 

In addition to the near and far 3-kpc arm population we find the maser l,v distribution demon- 
strates a prominent feature near —22° longitude which we interpret as the approximate tangent 
of the 3-kpc arms. Given the association of 6.7-GHz methanol masers with the parallel sections 
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of the 3-kpc arms and the tangent, we can expect the entire structure to exhibit methanol maser 
emission and form an approximately oval structure in the l,v domain. When this structure is com- 
bined with the spiral arms and other features, essentially all the maser emission is accounted for. 
This oval structure can either relate physically to an elliptical ring or a circular ring with radial 
motion ('expansion'). The average surface density of sources in the Galactic plane, assuming an 
annulus of 0.5 kpc centred on 3.5 kpc, is 6 and 7 sources per kpc 2 for the near and far portions of 
the arms (±15° longitude) respectively. This is comparable to the expected 'smeared' density of 
sources in the spiral arms ([Caswell et al.ll2010l ). The average density along the tangent (excluding 



sources associated with other structural features, i.e. limiting the velocity range to — 60kms 1 to 
-85kms _1 ) is 22 sources per kpc 2 . This higher density can be accounted for by the presence of a 
spiral arm origin, which we examine in Section | 



The existing models of the dynamics on the inner Galaxy account for <34% of the 3-kpc arm 
sources. In comparison, a simple elliptical ring or 'expanding' circular ring account for of the order 
of 95%. Although there are several models for an 'expandi ng' circular ring close to the centre of 



our Galaxy (with in a Galactocentric radius of <1 kpc) (e.g. IScoville 



1972 



Kaifu et al 



1996 



Ball y et al.lll987f). and nume rous examples of barred external galaxies with rings (e.g. 



197: 



2. 



Buta 



197 



4; 



198 



6. 



Treuthardt et al.ll2009l ). dynamical models of the Milky Way since the early examples in the 



1970s do not favour such a feature. Some radial motion may be a result of the influence of the 
radius of corotation resonance: a product of the transition between the elliptical orbits which follow 
the bar inside the radius of corotation resonance and the circular orbits which exist outside, with 
the 3-kpc ring the interface region betwee n the two, p ulling elliptically orbiting material into a 
circular orbit (an early example is that of IShandll972l ). The more favoured picture (in terms of 
dynamical modelling) of the elliptical ring or stream lines can account for the parallel sections and 
to some degree the negative tangent, but cannot readily account for the positive longitude masers 
with large velocities. The 6.7-GHz masers, as tracers of the structure, have the potential to fully 
define the spatial and kinematic structure of the 3-kpc ring, since further astrometric measurements 
will allow parallax distances (thus a direct measure of the ellipticity of the ring) and the full three 
dimensional velocity behaviour. 

The spiral arms are believed to originate at approximately the radius of cor otation resonance of 
the Galactic bar (e.g. Lopez-Corredoira et al. 1999 : Englmaier fc Gerhardlll99"9 ). where the Galactic 
material is rota ting at the same speed as the pattern speed of the bar, and radial oscillation is zero 
(|Lindbladlll974h , Highly elliptical orbits exist within the radius of corotation resonance, in contrast 
to circular orbits outside ( Englmaier k, Gerhardll2006 ). This requires that the 3-kpc arms should be 
either a mor e elliptical structure interior to the radius of corotation re sonance, with orbits following 
the bar (e.g. IContopoulos Papavannopoulos!ll980l : lEnglmaierl 12000) or a more circular structure 
close to the radius of corotation resonance, similar to resonance ring s seen in external galaxies, such 



as NGC2523 and NGC4245 (e.g. |ButaJ ll99J ; |Treuthardt et alJl200a ). As discussed in Sections E2H 



and 13.2.21 the maser distribution leads us to believe the long thin bar is the most influential on star 
formation, with a bar semi-major axis of ~3.4 kpc and orientation of ~45°, and the 3 kpc arms 
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appear to be well approximated by an elliptical ring (or circular with radial velocity components 
resulting from the resonance). This indicates a radius of corotation resonance of ~4.0 kpc, and a 
possible location for the origins of the spiral arms. 



3.3. The spiral arm origins 



Although th e tangents of all the spiral arms have been recognised in CO, Hi and other tracers 
such as CS (e.g. lBronfmanlll992l ) and IR star counts (e.g. iBeniamin et al.ll2005l ). identification of 
the arm origins, which fall within the terminal velocity envelope rather than at its boundary where 
the tangent points lie, has been elusive. This has meant models generally incorporate arbitrary 
or inferred arm origins within the l,v domain: the two 'major' arms (Perseus and Crux-Scutum) 
originate at opposite ends of the Galactic bar and the two 'minor' arms (Norma and Carina- 
Sagittarius) are anywhere between a few degrees and 90° offset, either starting from the radius 
of corotation resonance directly or branching from the two main arms. However, studies of other 
Galaxies and dynamical models show that the major arms do not necessa rily start directly at the 
bar end, but can be seen offset by several degrees of Galactic azimuth (e.g. ISandage &; Bedkdll994l : 
Russeillhooah . 



3.3.1. Spiral arm induced density enhancements 



As a product of interstellar material undergoing the influence of density waves, spiral arms are 
regions of dense material and star formation. External galaxies show the spiral arms start close 
to the end of the bar or a ring surrounding a bar and, as already mentioned, dynamical models 
place the arms starting at the radius of corotation resonance (~ 4 kpc). These regions are likely to 
experience significant compression of material such that star formation, induced by density changes 
in the ISM, is likely to be enhanced. The remaining high density regions of masers within ±28° may 
thus be explained by the origins of the spiral arms. The exception is between longitudes +6° to —6° 
where all the spiral arms pass and there are sources on the solar circle (zero line-of-sight velocities) 
together with sources associated with the Galactic Centre Zone. In this region an overabundance 
is to be expected even in the absence of spiral arm origins. 

High positive velocities at longitudes 25° to 35° have been considered the approximate start of 
the Cru x-Scutum arm, with a region of massive star formation inferred from a sequence of IR obser- 



2000; 



vatiqns (LH ammer slev et al 



Cabrera-Lavers et al 



Clark et al 



1994 



Garzon et al.lll997l : lL6pez-Corredoira et al . 1999; Hammerslev et al 



200 80 and a m assive young cluster of red supergiants (jDavies et al 



2007 



20oJ I g et al.ll2010l ). We see a high density of methanol masers in this region 

(1 bin in excess of 5a po i SSOn in FigureE]). The Crux-Scutum arm is the only arm with velocities 
above 60kms _1 at these longitudes, it is known to have a tangent at ~35°, and originate at longi- 
tudes smaller than this. The lower estimate of the radius of corotation resonance is 3.5 kpc, which 



Sun 



Sun 



Fig. 5. — Schematic of the two simple interpretations of the 3-kpc arms. The blue represents the 3-kpc 
arms as either an elliptical structure or a ring. The black line is the Sun-Galactic Centre. Black short dashes 
represent the long bar and grey long dashes represent the short bar. 



Table 2: Summary of comparison of maser distribution with models of the inner Galaxy. De- 
tails are given for the six models together with percentage association wi th the 3 kpc arm mase r 



distribution. 
E99 



CR is the corotation resonance. References are: M86 



Bissantz et al.l (|2003l ) ; R08 = iRodriguez-Fernandez Combes 



Muld er fc Lieml (jl986n 



200 



Enelmaier fc Gerhard! (11999T): F99 = iFuxi (11999T): S99 = ISevenster et al.1 (|1999^ : B03 
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33 
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3.4 


20 
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20 
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3.4 


16 


4 inner arms 


F99 


25 


3.5 


4.0-4.5 
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is tangential at 24°. Hence the arm theoretically originates between 24° and 35° and we interpret 
the dense cluster of masers at ~26° as observational evidence of the start of this arm. Assuming 
the longitudes 24° and 35° mark the boundary of this region, and that it is a portion of an annulus 
with a thickness of 1 kpc centred at the corotation resonance radius of 4 kpc, the origin has a 
density of 11 sources pe r kpc 2 . The discovery of the massive young red super giant cluster in th is 



region (jFiger et al.ll2006l ) implied substantial starburst activity 10-20 Myr ago (jDavies et al.ll2007l ). 
but the significant presence of 6.7-GHz methanol masers indicates there is still substantial ongoing 
star formation. 

As mentioned in Section [3.2.11 the methanol maser distribution demonstrates a dense ridge of 
sources near I = —22° which is coincident with the negative longitude tangent of a ~3.5 kpc ring. 
The density of masers within this region is 22 sources per kpc 2 (Section I3.2.2H compared to the 
average of both the 3 kpc ring and spiral arms which is 6 or 7 sources per kpc 2 . This leaves 15 or 
16 sources per kpc 2 which could be accounted for with the origin of a spiral arm. The Perseus arm 
is the only spiral arm in conventional models extendable to this region and is likely to correspond 
to the most negative velocities of the ridge. 

We now turn to the minor arms, Carina-Sagittarius and Norma. From fitting a logarithmic 
spiral to objects associated with the Carina-Sagittarius arm, and extrapolating towards the Galactic 
centre, the origin of this arm i s often locate d within the fourth Galactic quadrant within about 



10° of the Galactic centre (e.g. iRusseill 120031 ) . We observe a high density clump near longitude of 
—8° (1 bin above 5a po i sson in FigureE]). This is either a very dense and active region of high mass 
star formation within the Norma arm (at non-tangential longitudes and more negative velocities 
than expected) or it is more likely the indication of the start of the Carina-Sagittarius arm and 
the overlapping of the l,v parameters of this arm with those of the Norma arm (with the further 
implication that the Norma arm branches from the Perseus arm). The starting point of the Norma 
arm is less clear. It is often speculated to originate in the first Galactic quadrant, ap proximately 



opposite to the other minor spiral arm Carina-Sagittarius (e.g. lOrtiz h Lepind Il993l ). We see a 
dense cluster of sources from longitudes +10° to +12° extending to positive velocities beyond the 
Carina-Sagittarius and Crux-Scutum arms, which could be interpreted as the start of the Norma 
arm. If the Norma and Carina-Sagittarius arms branch from the major arms outside the radius of 
corotation resonance, rather than starting near the ring at the radius of corotation resonance, the 
conditions of the ISM may be less compressed and the level of star formation will not be as high 
as the origins of the major arms (the branching scenario is a more logical scenario geometrically if 
the ring is elliptical). 

The associations of the density enhancements are summarised in Table Q] together with an 
approximation of the density of sources per kpc 2 . We identify the high density regions of masers 
not associated with the bar or the 3-kpc arms as enhanced star formation at the starting points 
of the spiral arms: Crux-Scutum at approximately +26° (velocities circa 100 kms -1 ); Norma at 
approximately +12° (velocities circa 50kms _1 ); Carina-Sagittarius at approximately —8° (veloci- 
ties circa — 55kms _1 ); and Perseus at approximately —22° (velocities circa — 55kms _1 ). We find 
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the Crux-Scutum and Carina- Sagittarius arm origins to have densities of 10-15 sources per kpc , 
similar to the densities seen at the ends of the Galactic bar. The densities of the Norma and Perseus 
arms are confused by additional arms lying on the line-of-sight and the 3-kpc arm tangent. The 
masers contained within the spiral arm origins have a median latitude of 0.0±0.1°, are typically 
contained within a range of about 30kms -1 and have median fluxes in the range 6 to 10 Jy. 



4. Conclusion 



We have examined the distribution and density of 6.7-GHz methanol masers in the longitude- 
velocity domain. Both the longitude-velocity distribution and its structure function demonstrate 
the presence of structures on small (<0.03°) and large (>3°) scales. Through smoothing the 
density distribution on the small scales and binning on the large scales we identify statistically 
significant dense regions of masers indicative of enhanced high-mass star formation within Galactic 
scale structures. The maser distribution supports the presence of a long thin bar inclined at angle 
of ~45° to the Sun-Galactic centre line of sight. The lack of methanol emission from a short, 
bulge associated bar implies the two Galactic bars represent different ages: the short boxy bulge 
bar an older structure and the long thin bar a younger structure, undergoing current high-mass 
star formation. We identify a prominent tangent of the 3-kpc arms near —22° and find the maser 
distribution of the 3-kpc arms is readily associated with a continuous ring structure in the longitude- 
velocity domain. High densities of masers identify the approximate starting points of the spiral 
arms: the major arms, Crux-Scutum at +26° and Perseus at -22°, slightly offset from the bar 
ends; the minor arms, Norma at 12° and Carina-Sagittarius at -8°, possibly branching from the 
major arms. The 3-kpc arm ring and spiral arm origins, combined with the spiral arms themselves, 
account for essentially all the 6.7-GHz methanol masers and their density enhancements within 
longitudes ±28°. 6.7-GHz methanol masers clearly delineate many of the important structures of 
our Galaxy and provide a new observational basis to constrain dynamical models. Multi- wavelength 
studies and astrometric distances through VLBI will further enhance the importance of this species 
of maser as a tool in understanding the structure of our Galaxy. 

Facilities: Parkes (Methanol Multibeam). 



A. Minor implications of uncertainties in spiral arm shape and rotation parameters 



As noted in Section [2] we adopt the lTaylor Corded (|1993l ) spiral arms and the lBrand 8z Blitz 



(|l993l ) rotation curve, together with a Galactocentric solar distance (Rq) of 8.4 kpc and a circular 
rotation of the Sun (Qq) of 246kms _1 . We now discuss the potential implications of uncertainties 
in our assumptions regarding the shape of the spiral arms and the rotational parameters of the 
Galaxy. 
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Fig. 6. — The elliptical stre amlines corresponding t o the 3-kpc arm (top) from IPeterd (jl975T ) and the 
expanding ring (bottom) from Cohen fc Daviesl ( 1976I) . The top is based on two streamlines, one with a 
semi major axis of 4.9 kpc with an ellipticity of 0.69 and a fixed angular momentum of 770 km s -1 kpc -1 , 
the other with a semi major axis of 4.6 kpc with an elliptiicity of 0.68 and a fixed angular momentum of 
680 km s" 1 kpc -1 . The bottom is based on a radius of 4.0 kpc and a rotational velocity of 210 km s -1 with 
an expanding velocity of 53kms _1 . 
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Fig. 7. — Fraction of maser sources consistent with the 3-kpc ring model for different parameters of the 
ellipse. Colour scale represents 10 to 100% source association. Left figures are for the association of 3-kpc 
arm sources only, right figures are for association with both the 3-kpc arm sources and the unassociated 
sources described in the text (with equal weighting to both). The top figures are for fixed semi major and 
semi minor axes, 4.1 kpc and 2.214 kpc respectively. The middle figures are for fixed semi-minor axis of 
2.214 kpc and fixed angular momentum of 320 kms -1 kpc -1 . The bottom figures are for a fixed orientation 
of 38° and angular momentum of 320 kms -1 kpc -1 . 
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Fig. 8. — Fraction of maser sources consistent with the 3-kpc ring model for different parameters of the 
ring. Colour scale represents 10 to 100% source association. Left figures are for the association of 3-kpc arm 
sources only, right figures are for association with both the 3-kpc arm sources and the unassociated sources 
described in the text (with equal weighting to both). The top figures are for a fixed outward Galactoccntric 
radial velocity of 48kms _1 , the middle figures are for a fixed circular velocity of 201kms _1 and the bottom 
figures are for a fixed radius of 3.4 kpc. Black lines delineate likely boundaries for parameters from models 
and observations (see text for details). 
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Fig. 9. — Longitude- velocity plot with a 3 kpc ellipse and ring. Symbols are the same as Figure Q] The 
ellipse is 0.5 kpc thick in the radial direction and has a semi-major axis of 4.1 kpc, semi minor axis of 2.2 
kpc, orientation of 38° and an angular momentum of 320 km s -1 kpc -1 . The ring is 0.5 kpc thick in the 
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Fig. 10. — Comparison of the 6.7-GHz methanol masers with model inner Galaxy structures. Symbols are 
the same as Figure CD Th e blue l ines delineate the regio n identified in CO emission as the 3-kpc arms by 
Dame fc Thaddeusl |2008l) . TOP: Sevenster et al. ( 1999) elli ptical ring with a bar orientation of 44° and a 
corotation radius of 4.5 kpc. MIDDLE: iBissantz et al.l ([20031 ) combi nation of elliptically orbiting gas with a 
bar orientation of 20° and a corotation radius of 3.4 kpc. BOTTOM: Rodriguez-Fernandez fc Combes ( 20081) 
lateral arms with a combination bar of orientations of 20° and 75° and a corotation radius of 5-7 kpc. 
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Fig. 11. — Comparison of the 6.7-GHz methanol masers with model inner Galaxy structures continued. 
Symbols are t he same as Figure [0 The blue li nes delineate the region identified in CO emission as the 
3-kpc arms by Dame fc Thaddeus ( 20081 ). TOP: Mulder fe Liem ( 19861 ) four inner arm-like features: near 
and far inner lindblad radius arms and i nner ultra-harmonic arms. T his model has a 20° bar orientation 
and corotation radius of 4kpc. MIDDLE: Englmaier fc Gerhard ( 1999h four inner arm -like featur es from the 
bar ends with a bar orientation of 20° and a corotation radius of 3.4 kpc. BOTTOM: iFuxl (|1999l ) four inner 
arm-like features: the 3-kpc arm, the 135 km s -1 arm and the near and far connecting arms. This model 
has a bar orientation of 25° and a corotation radius of 4— 4. 5 kpc. 
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A.l. Choice of spiral arm shape 



Our knowledge of the Galactic location and spatial properties of the spiral arms has been 
limited by the accuracy of the distances of objects tracing the spiral arms , such as Hn regions and 
Giant Molecular Clouds. The often cited model of lGeorgelin Georgelinl (|1976l ) is based on inter- 
pretation of star forming regions with assigned distances, many of which are kinematic distances, 
known to suffer significant error margins (see for example iGomea 120061 : iReid et al.1 120091 ) . Models 
are often made to fit the observed tangential velocities of atomic (Hi) and molecular (CO) emission, 
interpolating the structure between these extremes to approximate what is seen observationa l ly (e.g . 
Baba et al.ll2010l ). As such, widely used models of the spiral arms, such a s Tavlor fc Corded (119931 ). 
have changed minimally in the past 44 years since Georgelin Sz Georgelin ( 197(jI ) , so suffer the same 
uncertainties and some caution should therefore be taken with the location of the spiral arm loci. 
Statistical arguments summarising the mod els suggest there is a fay o ured picture of a logarithmic 
four arm spira l wit h a pitch angle o f ~12° |Valleelll995l . I2OO2I . bood . hoosh . The recent studies of 
Russeill (|2003l ) and iHou et al.l (|2009j), with logarithmic and polynomial fits respecti yely, favour a 



four a rm model and specifically one which has two 'major' and two 'minor' arms. iLevine et al 



(|2006l ) found concurring evidence for a four arm structure in the Hi emission in the outer Galaxy. 
Additionally, although the Galaxy appears to have an overall grand structure, spurs, such as the 
local 'Orion arm', exist. 

It is the combination of observations at multiple wayelengt hs and the 'piggybacking' of dis- 
tances through associations of these (e.g. lUrquhart et al.ll2010i ). that is helping to improve our 
current understanding and models. Ultimately however, the thousands of stellar distances of the 
upcoming NASA Gaia mission together with astrometric parallaxes of masers through VLBI (iden- 
tifying the optically obscured star forming regions) will have the potential to truly delineate the 
spiral arms. Unfortunately, both of these are several years away, making the loci of models our 
current best estimate: until then the likely overall structure, number of arms and continuity of the 
arm loci will remain helpful in recognising potential spiral arm association of maser sources. 



A. 2. Choice of rotation curve 



As the lTavlor h Corded (|1993l ) arms are largely based on those of lGeorgelin h Georgelinl (|1976l ) 
and these were based on kinematic distances using a rotation curve, it is pertinent to adopt a similar 
rotation curve to convert the arm shapes to the l,v domain. iGeorgelin &; Georgelinl (| 19761 ) used a 
rotation curve with a Q of ~250 km s -1 and an Rfp of 10 k pc. When scaled to 8.5 kpc this equates to 
a O© of ~220kms~ 1 , for which the lBrand Sz Blita ( 19931 ) rotation curve is a goo d approximation . 
However, we also tested the association of ma s ers w ith th e rotation c urves of IClemend (j 19851 ) , 
Fich et al.l (|l989h . iMcClure-Oriffiths fc Picked (Eoo3) and faeid et al.l (hood) . With lkpc thick 



spiral arms and a 7kms 1 velocity tolera nce in assoc i ation we find: 



77% of the sources with the spiral arms; iFich et al.l (|1989l ) linear curve 78%; iFich et al.l (jl989j) 



Clemend (1985) accounts for 
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power -law curve 79%; iMcClure-Griffiths & Dickevl (|2007l ) 80%; and the flat rotation of iReid et al 



(|2009l ) accounts for 76% of the sources. Although individual examples of associated sources change 
between the rotation curves, the total remains in the range 76—80%, mirroring the conclusions of 
the previous section and emphasising that the inferred total number of spiral arm sources is largely 
unaffected by rotation curve choice. 



A. 2.1. Solar Parameters 



In addition to the choice of rotation curve, Galactic kinematics are dependent on the parameters 
of solar motion and the local standard of rest. For the solar distan ce, the con sensus appears 



to be 8.4 kpc, based on or bits of SO-2 stars implying 8.4±0.4 kpc (jGhez et all 120081 ') and 8.33±0.35 
kpc (jGillessen et al.ll2009l ). T his is consistent w ith the parallax observations of water masers in 
Sgr B2 indic ating 7.8±0 . 8 kpc (jReid et al.ll2009l ) and the parallaxes of methanol masers indicating 
8.4±0.6 kpc (jReid et alibood ). 



The IAU standards of solar motion are implicitly incorporated in our maser LSR velocities. 
However, recent maser astrometric observations have prompted revision. Two of the parameters, 
solar motion towards the Galactic centre, U®, and solar motion towards the north Galactic pole, 
Wq, have remained largely unchanged (10.3 kms -1 and 7.7 kms -1 respectively), but the solar mo- 
tion in the direction of Galactic rotation, Vq, has not. Originally Vq was chosen to be 15.3 kms -1 
(the IAU standard), then it was revised to 5 .25±0.62 kms" 1 based on stellar kinematics of the 
Hipparcos catalogue (jDehnen Binnevl I1998TI. However the d ataset was re-examined, resulting 
in further revision to 12.24±0.47 kms -1 ( Schonrich et al,ll2010l ). This higher value of Vj^, return 



i ng to a v alue close to the orig i nal IA U standard, is in line with the suggestions of IReid et al 



(|2009l ) and McMillan Binneyi (120101 ) based on maser proper motions. The revision to a higher 
value may account for th e apparent motion of star forming regions counter to Galactic rotation 
(|Bobvlev fc BaikovJl2010h . 



Galactic circular rotation of the Sun, Qq has also been subject to recent revision. The IAU 
standard of 0©, 220 km s -1 , has been revised to a higher value based on the results of proper motions 
de duced from maser astrometry observations. The revised values are 254zfcl6k ms~ 1 as e stimated 



by lReid et all (j2009h . or either 246±30kms -1 or 244±13kms -1 as estimated by lBovy et all (|2009l ) 



Therefore the current best estimates of the parameters are: Uq = 11.1 kms -1 ; Wq = 7.25 kms -1 
Vq = 12.2 kms -1 ; Qq = 246 kms -1 ; and Rq = 8.4 kpc. The most significant difference is the in- 
crease in ©q. We adopt the Rq of 8.4 kpc throughout, but explored the impact of the other 
variations on our conclusions, and the effect is minimal. The adjustment of @q has a reasonable 
impact on kinematic distance estimates, but as we are working in the l,v domain, it does not in- 
fluence our maser parameters. It does affect the location of the spiral arm loci if we transfer them 
with this value rather than the IAU standard, but this only serves to shift the loci in velocity by 
~ 5-10 kms - . 
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